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Abstract. We investigate a class of generalized Schrodinger operators 
in L 2 (R 3 ) with a singular interaction supported by a smooth curve 
r. We find a strong-coupling asymptotic expansion of the discrete 
spectrum in case when T is a loop or an infinite bent curve which 
is asymptotically straight. It is given in terms of an auxiliary one- 
dimensional Schrodinger operator with a potential determined by the 
curvature of T. In the same way we obtain an asymptotics of spectral 
bands for a periodic curve. In particular, the spectrum is shown to 
have open gaps in this case if T is not a straight line and the singular 
interaction is strong enough. 

1 Introduction 

The subject of this paper are asymptotic spectral properties for several classes 
of generalized Schrodinger operators in L 2 (K 3 ) with an attractive singular 
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interaction supported by a smooth curve or a family of such curves. On a 
formal level, we can write such a Hamiltonian as 



however, a proper way to define the operator corresponding to the formal 
expression is involved and will be explained in Sec. 12.21 below 1 . A physical 
motivation for this model is to understand the electron behavior in "leaky" 
quantum wires, i.e. a model of these semiconductor structures which is re- 
alistic in the sense that it takes into account the fact that the electron as a 
quantum particle capable of tunelling can be found outside the wire - cf. |EI 
for a more detailed discussion. 

One natural question is whether in case of a strong transverse coupling 
properties of such a "leaky" wire will approach those of an ideal wire of zero 
thickness, i.e. the model in which the particle is confined to T alone, and how 
the geometry of the configuration manifold will be manifested at that. In 
the two-dimensional case when T is a planar curve this problem was analyzed 
in |KY1| IEY2j where it was shown that apart of the divergent term which 
describes the energy of coupling to the curve, the spectrum coincides asymp- 
totically with that of an auxiliary one-dimensional Schrodinger operator with 
a curvature-induced potential 2 . 

The case of a curve in M 3 which we are going to discuss here is more 
complicated for several reasons. First of all, the codimension of T is two 
in this situation which means that to define the Hamiltonian we cannot 
use the natural quadratic form and have to employ generalized boundary 
conditions instead. Furthermore, while the strategy of |EY1| IEY2j based on 
bracketing bounds combined with the use of suitable curvilinear coordinates 
in the vicinity of T can be applied again, the "straightening" transformation 
we have to employ is more involved here. Also the bound on the transverse 
part of the estimating operators are less elementary in this case. 

Let us review briefly the contents of the paper. We begin by constructing 
a self-adjoint operator H ar which corresponds to the formal expression (jl.ljl . 
where T is a curve in M. 3 ; this will be done in Sec. 12.51 To this aim we 
employ in the transverse plane to T the usual boundary conditions defining 

1 In particular, this is the reason why we use here a formal coupling constant different 
from the parameter a introduced in the condition l|2.4|l below. 

2 A similar analysis was performed in |Exj for smooth surfaces in R 3 where the asymp- 
totic form of the spectrum is given by a suitable "two-dimensional" operator supported 
by the surface T. 



A - a5(x - T) , 
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a two-dimensional point interaction [AGHH, Sec. 1.5]. Recall that the latter 
is known to have for any a G K a single negative eigenvalue which equals 
£ a = _4e 2 (- 2 ™+^( 1 )), where -^(1) = 0.5777... is the Euler constant. The 
main topic of this paper are spectral properties of H a ^ in the strong-coupling 
asymptotic regime which means here that —a is large. The auxiliary operator 
mentioned above is given by 

S: =-A--K 2 , 
4 

where A is the one- dimensional Laplace operator on the segment parame- 
terizing T and k is the curvature of T. Its discrete spectrum is non-empty 
unless T is a straight line; we denote the j-th eigenvalue as fij. Our main 
results can be then characterized briefly as follows: 

Discrete spectrum: If T is a loop, we show in Sec.EJthat the j-th eigenvalue 
Aj(a) of H a ^r admits an asymptotic expansion of the following form, 

Xj(a) = £ a + fij + 0(e na ) as a — ► — oo , 

and the counting function a i— > #<Td(if Qi r) satisfies in this limit the relation 

#a d (H a ) = -(-Q^(l + 0(e™))- 

7T 

In addition, the last formula does not require T to be a closed curve as 
we shall show in Sec. 13.51 Moreover, if T is infinite with k ^ and at the 
same time asymptotically straight in an appropriate sense then the above 
expansion for \j(a) holds again - cf. Sec. |U 

Periodic curves are discussed in Sec. El we perform Bloch decomposition 
and use the same technique as above to estimate the discrete spectrum of 
the fiber operators. In particular, we find that if T is periodic curve and k(-) 
is nonconstant then a(H a -p) contains open gaps for —a sufficiently large. In 
the closing section we will show that the problem can be rephrased in terms 
of a semiclassical approximation and list some open problems. 
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2 Hamiltonians with curve-supported 
perturbations 

2.1 The curve geometry 

Let T be a curve in IR 3 (either infinite or a closed loop) which is assumed to be 
C k , k > 4. Without loss of generality we may assume that it is parameterized 
by its arc length, i.e. to identify F with the graph of a function 7 : / — > M. 3 , 
where / = [0, L] (with the periodic boundary conditions, 7(0) = •y(L) and 
the same for the derivatives) if T is finite and I = M otherwise. One of our 
tools will be a parametrization of some neighbourhoods of T. To describe it 
let us suppose first that the curve possesses the global Frenet's frame, i.e. 
the triple (t(s),b(s),n(s)) of tangent, binormal, and normal vectors which 
are by assumption C k ~ 2 smooth functions of s G /; recall that this is true if 
the second derivative of T vanishes nowhere. 

The mentioned neighbourhoods are open tubes of a fixed radius centred 
at T: given d > we call fl d := {x G M 3 : dist(x,T) < d}. We will 
impose another restriction on the class of curves excluding those with self- 
intersections and "near-intersections", i.e. we suppose that 

(an) there exists d > such that the tube fl d does not intersect itself. 

Our aim is to describe fid by means of curvilinear coordinates, i.e. to write it 
as the image of a straight cylinder Bd '■— {r G [0, d), 6 e [0, 2n)} by a suitable 
map. If the Frenet frame exists we choose the latter as (pd '■ Dd — ► K 3 defined 



Ms, r, 0) = j(s) - r [n{s) cos(0-/3(s)) + b(s) sin(0-/3(s))] , (2.1) 



where Vd := I x Bd and the function (3 will be specified further. For conve- 
nience we will denote the curvilinear coordinates (s, r, 9) also as q with the 
coordinate indices (1, 2, 3) <-> (s, r, 8), and moreover, since it can hardly lead 
to a confusion we use the same notation (pd for the mappings with target 
spaces M 3 and fid which will need later. 

The geometry of fid is naturally described in terms of its metric tensor 
(gij); the latter is according to |DEj expressed by means of the curvature k 
and torsion r of T in the following way 



by 




4 



where 



<; := r — /3 S and h := 1 + rncos(9 — (3) . 



(2.2) 



We use here the standard conventions (3 >s = <9 s /3 and gr y ' = (<7ij) _1 . In partic- 
ular, the volume element of is given by <if2 = g x l 2 dq where g := det(^jj). 
The simplest situation occurs if we choose 



because then the tensor takes the diagonal form = diag(/i 2 , l,r 2 ). 

Remarks 2.1 (a) It is well known that compact manifolds in M n have the 
tubular neighbourhood property. Thus if T is a finite C 4 curve then the 
assumption (an) is satisfied iff T has no self-intersections. 

(b) Combining the explicit formula for g^ with the inverse function theorem 
it is easy to see that the inequality d||/c||oo < 1 is sufficient for (p d to be locally 
diffeomorphic. 

The special rotating system described above is called in the theory of 
waveguides usually Tang system of coordinates. If I is finite, the functions 
h, h, h are bounded by assumption, while in the case I — M. the global 
boundedness has to be assumed. The main problem, however, is that the 
described construction may fail if the Frenet frame is not uniquely defined. 
Hence we suppose in general that 

(ar2) for all d > small enough there is a diffeomorphism 0^ : T>d — > Qd 
such that the corresponding metric tensor is = diag(/i 2 , 1, r 2 ) where 
h is given by (|2.2|) with j3 which is locally bounded, C k ~ 2 smooth with 
a possible exception of a nowhere dense subset of J, and h together 
with its first two derivatives are bounded. 

While it represents a nontrivial restriction, this hypothesis can be neverthe- 
less satisfied for a wide class of curves without a global Frenet frame. 

Example 2.2 Suppose that the curve parameter interval / can be covered 
by at most countable union Uj g j Ij °f intervals Ij = [aj, bj] such that a pair of 
different Ij, Ik has in common at most one endpoint, and furthermore, either 
the Frenet frame exists in (aj,bj) or T^ ss = in [a,j,bj]. In the former case 
we assume also that limits of n(s) as s approaches a,j and bj exist. We claim 
that in such diffeomorphism with a diagonal can be constructed. 



(2.3) 
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Let us describe first its building blocks. If the Frenet frame exist in (aj, bj) 
we construct a map (pf on the appropriate part of the curve by (|2.1j) with (3 
replaced by a function (3j satisfying the condition ()2.3|) . On the other hand, 
in the straight parts we construct <j>^' similarly choosing a constant for f3j 
and an arbitrary pair of unit vectors forming an orthogonal system with the 
tangent for b, n with a proper orientation. 

The maps <f>^ can be patched into a global diffeomorphism by choosing 
properly the /3/s. Suppose first that the family {cij} of left endpoints has 
no accumulation points in the interior of I. In that case we may identify 
without loss of generality the index set J with a segment of Z, i.e. j = 
-N, -N + 1, . . . , M for some N, M e N U {00} and to suppose that the 
interval family is ordered, a n+ \ = b n . Assume now that such a diffeomorphic 
map exists on \J™ = _ n Ij- The left and right limits of the base vector systems 
at the points a_ n = 6_ n _i and b m = a m+ i exist by assumption and differ at 
most by a rotation in the normal planes to T at these points, so the map 
can be extended to a diffeomorpism on Uj^t-n-i Ij by adjusting /3_ n _i and 
(3 m+ i, notice that the condition ()2.3|) remains valid on (aj, bj) if the function 
(3j is shifted by a constant. The sought conclusion then follows easily by 
induction. On the other hand, let the set {aj} have accumulation points in 
the interior of I; we call them {c^} ordering them into a increasing sequence, 
finite or countable. The above construction defines a diffeomorphism in each 
interval Ik := (ct, Ck+i) between adjacent points. Then the argument can be 
repeated, the role of Ij being now played by the intervals Ik- 

Having constructed such a <pd : T>& — > fid one can check directly whether 
the global boundedness conditions of the assumption (ar2) are satisfied. 

2.2 Singularly perturbed Schrodinger operators 

The Hamiltonians we want to study are Schrodinger operators with s-inde- 
pendent perturbations supported by the curve T. Such operators can be 
understood as the Laplacian with specific boundary conditions on T and the 
aim of this section is to make this conditions precise. 

Let us assume that for a number d the map (f>d satisfies conditions (an, 2). 
Given p 6 (0, d) and #0 G [0, 27r) denote by T p ^ the "shifted" curve lo- 
cated at the distance p from V which is defined as the (pa image of the set 
I x {p, # } c V d ; recall that the global diffeomorhism <p d exists by assump- 
tion (ar2). Consider the Sobolev space W£T(A \ T), where A is an open 
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bounded or unbounded set in R 3 such that Qd C A; particularly A may co- 
incide with whole R 3 . Since its elements are continuous on A away of T, the 
restriction of a function / £ W^(A \ T) to the "shifted" curve located suffi- 
ciently close to T is well defined; we will denote it as / lr p So (')• I* 1 f ac ^ we can 
regard / \ r g as a distribution from D'(0, L) parameterized by the distance 

p and the angle 6 . We shall say that a function / £ W^(A \ V) n L 2 (A) 
belongs to T^ d if the following limits, 

Q(f)(s) := lim \f\ rpeo (s) + E(f)(s)lnp 



exist a.e. in [0, L], are independent of 9 , and define a pair of functions 
belonging to L 2 (0, L); for an infinite curve [0, L] is replaced by R. We should 
also stress here that the elements of (A\T) are in fact distributions from 
^'(R 3 ), however, in the definition of Ta we can naturally identify them with 
their canonical imbeddings into L 2 (A). 

Given a function / £ Ta we write f~a.bc(T) if the limits S(/)(-), fl(f)(-), 
characterizing the behavior of / close to T satisfy the following relation 

27raS(/)( a ) = n(/)(a). (2.4) 

With these prerequisites we can define the singularly perturbed Schrodinger 
operator in question through the set 

D(H a>v ) = {/ £ T R3 : f~a.bc{T)} 

on which the operator H a p : D(H a ^) — > L 2 (R 3 ) acts as 

H a ,rf(x) — — A/(x) , x£R 3 \T. (2.5) 

To show that if Qj r makes sense as a quantum mechanical Hamiltonian we 
will assume here that T is finite or infinite periodic. Another interesting case, 
that of an infinite non-periodic curve which is asymptotically straight, needs 
additional assumptions and will be discussed separately in Sec. HI 

Theorem 2.3 Under the stated assumptions H a r is self-adjoint. 
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Proof. One check using integration by parts and passing to the curvilinear 
system of coordinates q = (s, r, 9) in a sufficiently small tubular neighbour- 
hood of T that the following boundary form, 

v ■ v(f,g) = (H atr f,g) - (f,H a>r g) 

vanishes for all /, j £ D{H a ^) 1 i.e. that the operator H a r is symmetric. 
To check its self-adjointness we can proceed in analogy with |EK| Thm. 4.1]. 
Repeating the argument presented there step by step we derive the resolvent 
of H a Y and the sought result then follows from \Po\ Theorem 2.1]. An 
alternative way is to note that H a Y is one of the self-adjoint extensions 
discussed in |Kuj . It is true that in this paper stronger smoothness conditions 
for T were adopted, however, the results remain valid for the C 4 class. ■ 

The operator H a ^ will be a central object of our interest. It is natural to 
regard it as a Schrodinger operator with the singular perturbation supported 
by the curve T. 

Remarks 2.4 (a) The choice of boundary conditions (J2.4|) which we used in 
the construction had a natural motivation. If T is a line in M 3 one can sepa- 
rate variables; in the cross plane we then have the two dimensional Laplace 
operator with a single-centre point interaction — A Qj {o} which is a well stud- 
ied object - cf. AGHH, Sec. 1.5]. To define it, one considers for a function 
/ E W^(R 2 \ {0}) n L 2 (R 2 ) the following limits 

^ (/) := ~ H rnV ' ^ (/) := + m Inr) ; 
if they are finite and satisfy the relation 

27TqE(/) = Q(f) , (2.6) 

the function / belongs to the domain of -A Qi {o}. Using the explicit form of 
its resolvent it is easy to see that such an operator has for any a G M. exactly 
one negative eigenvalue which is given by 

f Q = _ 4e 2(-2™ + ^(i)) ; ^ = _ .577... (2.7) 

Obviously, it coincides with the bottom of the essential spectrum of H a>r 
for a straight T. We know from EKJ that this property is preserved if T 
is curved but asymptotically straight in a suitable sense; in that case the 
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operator has a non-empty discrete spectrum - cf. Sec. EJ It is also clear from 
the relation (|2.7|) and the corresponding eigenfunction [AGHH1 Sec. 1.5] that 
a strong coupling corresponds to large negative values of a. 

(b) For the sake of brevity we use in analogy with ()2.4|) for the boundary 
conditions (|2.6|) the abbreviation f~a.bc(0), later we employ similar self- 
explanatory symbols for other conditions, Dirichlet, Neumann, periodic, etc. 



3 Strong coupling asymptotics for a loop 



In this section we will discuss in detail the strong-coupling asymptotic be- 
havior of the discrete spectrum in the simplest case when T is a finite closed 
curve satisfying the regularity assumptions stated above; by Remark 12. lT a) 
it means that T is C 4 and does not intersect itself. 

Remark 3.1 In the following considerations we will rely on the operator 
inequality A < B, where both of operators A, B are self-adjoint and bounded 
from below. To be precise we are going to follow the definition from [RS, 
Sec. XIII. 15], i.e. A < B iff 

q A [f}<qB[f], f e Q{B) C Q(A) , 

where q^, Qb are the forms associated with A, B having the form domains 
Q(A),Q(B), respectively. 

Since T is compact it does not influence the essential spectrum of H a Y- 
This can be seen by writing explicitly the resolvent |Poj and checking that it 
differs fr om the free one by a compact operator in analogy with the argument 
used in BKKS for codimT = 1. However, there is a simpler way. 



Proposition 3.2 With the stated assumptions we have 

o- ess (H a: r) = a ess (-A) = [0, oo) . 

Proof. By Neumann bracketing we can check that inf o~ css {Ha,r) = 0. Indeed, 
choose a ball B such that T is contained in its interior and call H^f? the 
Laplace operator in L 2 (IR 3 ) with the same boundary condition on T as H a r 
and Neumann condition at dB. We have H a ^ > H^ B and the spectrum 
of the latter is the union of the interior and the exterior component. The 
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first named one is discrete and the spectrum of the other is the non-negative 
halfline, so the claim follows from the minimax principle. To show that every 
positive number belongs to a(H a ^) it is sufficient to construct a suitable Weyl 
sequence; one can use a Weyl sequence for —A chosen in such a way that its 
elements have supports disjoint from B. ■ 

Let us turn to the main subject of this section. To describe how the 
discrete spectrum of H a behaves asymptotically for a — > — oo we employ the 
comparison operator defined by 

S = -^-^-: D(S)^L 2 (0,L), (3.1) 

with the domain D(S) = {4> G W 2 ' 2 (0, L); 4>~p.bc(0, L)}, i.e. determined 
by periodic boundary conditions, 0(0) = 0(L),0'(O) = 4>'(L). Furthermore, 
«(•) is the curvature of T. It is worth to stress that S acts in a different 
Hilbert space than H a< r- We denote by fij the j-th eigenvalue of S. With 
this notations our main result looks as follows: 

Theorem 3.3 (a) To any fixed neN there exists an a(n) G R such that 

#cr d (H aj r) >n for a< a{n) . 

The j-th eigenvalue Xj(a) of H a T admits an asymptotic expansion of the 
following form, 

Xj(a) = (, a + fij + 0(e Ka ) as a — >• — oo . 
(b) The counting function a i— > #o~d{H a p) behaves asymptotically as 

# ( T d (^) = -(-ea) 1/2 (l + 0(e™)). 
7T 

The proof of the theorem is divided into several steps which we will describe 
subsequently in the following sections. It is also worth to stress here that the 
error term 0(e na ) is not uniform with respect to j; this will be clear from 
Lemma 13.51 below. 
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3.1 Dirichlet-Neumann bracketing 

Our aim is to estimate the operator H a p in the negative part of its spectrum 
from both sides by means of suitable operators acting in a tubular neigh- 
bourhood Qd of T, with d sufficiently small to make the assumptions (an, 2) 
satisfied. The first step in obtaining the estimating operators is to impose 
additional Dirichlet and Neumann condition at the boundary of f^. Let thus 
the operators H 3 aVl j = D,N, in L 2 (fl d ) act as the Laplacian with the do- 
mains given respectively by D(H 3 aT ) = {/ <E T Qd : f~a.bc(T), f~j.bc(dtt d )}; 

it is straightforward to check that operators H 3 aT are self-adjoint. Now the 
well-known result |RS| Sec. XIII. 15] says that 

- A% d © Hl T < H a , T < -A° d © H° T , := R 3 \Q d . 

What is important is that the operators — corresponding to the exterior 
of fl d do not contribute to the negative part of the spectrum because they 
are both positive by definition. 

It is convenient to express the operators H 3 aT in the curvilinear coordi- 
nates q = (s, r, 8); this can be done by means of the unitary transformation 

Uf = fofa: L 2 (Q d ) - L 2 (V d , g l ' 2 dq) , V d = [0, L] x B d ; 

recall that the global diffeomorhism <pd exists by assumption (ar2). Then 
the operators H 3 aT := UH 3 aT U~ l act as 

f(x)^-(g- 1 / 2 d l g 1 / 2 g l3 d 1 f)(x) for xen d \T 

with the domains {/ G T Qd : f~a.bc(T), f~j.bc(u r (d)), f~p.bc(u s (0),u s (L))}, 
respectively, where we have introduced the notation 

:= {q^V d :qi = t}. 

To simplify it further we remove the weight g 1 ! 2 appearing in the inner prod- 
uct of the space L 2 (V d , g x l 2 dq). This is done by means of the another unitary 
map, 

U : L 2 (V d , gV*dq) - L 2 (V d , dq) , Uf := g^f ; 

the images of H 3 aV will be denoted as H 3 aV = UH 3 aT U~ l . The aim of these 
unitary transformations is to find a representation where the eigenvalues - 
which we need to estimate the eigenvalues of H a Y by means of the minimax 
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principle - are easy to analyze. A straightforward calculation analogous to 
that performed in [DEj yields explicit formulae for H^ v , j = D,N, which 
both act as 3 

-d^dj - \r- 2 + V , 
where V is the effective potential given by 

V = g- 1/4 (d l9 v(d^)) + \r~\ (3.2) 
while their domains are different, 

D{HZv) = {f e?v d : g- 1//L f~a.bc(r), f~P-bcM0),u> s (L)), 

f~D.bc("r(d))}, 
D(H£r) = {feT Vd : g-Wf'a.bc(F), f~p.bcM0),u a (L)) , 

(d r f) r=d = -[{g^drg-^fUd} , 

Remark 3.4 Notice that the boundary conditions satisfied by functions 
from D(H J aT ) on the curve V can be written in a simpler way. Since only 
the leading term in g™ 1//4 is important as r — * 0, they are equivalent to 
r -1 / 2 f~a.bc(T). Notice also that while the Dirichlet boundary condition at 
dVtd persists at the unitary transformation, the Neumann one is changed by 
U into a mixed boundary condition. 

3.2 Estimates by operators with separated variables 

While the operators H 3 a r , j — D, N, give the two-sided bounds for the neg- 
ative eigenvalues of H a , they are not easy to handle. This is why we pass to 
a cruder, but still sufficient estimate by operators with separated variables. 

In the first step we will make the boundary conditions in the lower bound 
independent of the coordinates. The boundary term involved in the definition 
of D(H^ r ) depends on s and 9. We replace the corresponding coefficient by 
M := ||5 ,1 / 4 <9 r .g _1 / 4 || TOCI passing thus to the operator 

II 1 1 L (Lu r (a)) 

Kv ~ A ^ ® 1 + 1 ® (~ A a ) + v <Kr 

3 We employ the usual convention that summation is performed over repeated indices 
keeping in mind that (g lJ ) is diagonal. 
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on L 2 (0,L) ® L 2 {B d ), where -A h := -d s h~ 2 d s : £>(£) -> L 2 (0,L) and 

"A; := -<9 r 2 - r- 2 «9 e 2 - V 2 : D(A^ L 2 (5 d ) , 

D(A-) := {/G<' c 2 (5 rf \{0}): E L 2 (B d ), r^O), 

(<9 r /)| r=d = M/| r= 4 

with the boundary condition at the centre of the circle written in the simpli- 
fied form mentioned in Remark 13.41 The upper bound contains no boundary 
term depending on s or 9 so we can put 

Kr = H° r = -A h ®I + I® (-A+) + V 

which acts in the same way but the above mixed boundary condition on dB d 
is replaced by the Dirichlet condition. 

The next estimate concerns the effective potential V given by ()3.2|) ; by a 
straightforward calculation |DE| we can express it in terms of the curvature 
together with the function h and its two first derivatives with respect to the 
variable s as follows, 

It is important that up to an 0(d) term this expression coincides with the 
potential involved in the comparison operator S. Indeed, since h is continuous 
on a compact set and thus bounded, by (|2.2j) there exists a positive Ch such 
that the inequalities 

C£(d) < h- 2 < C+(d) with (d) := 1 ± C h d, 

hold for all d small enough. Since T is C 4 by assumption, the derivatives h iS 
and h )SS are also bounded; hence (|3.3|) yields the estimate 



with a positive Cy valid on for all sufficiently small d. At the same time, 
we can apply the above bounds for h~ 2 to the longitudinal part of the kinetic 
term. Putting all this together we get 

L7, ®I < -A h ®I + V<L+®I, 
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where 

^ := -Cf ^ - * ± C v d : £>(£) -> L 2 (0, L) . 

Summarizing the above discussion, we can introduce a pair of operators with 
the longitudinal and transverse components separated, namely 

B± :=L±®I + I® (-A±) on L 2 (0, L) g) L 2 (B d ) , (3.4) 

which give the sought two-sided bounds, ±-ff^ r < 



3.3 Component eigenvalues estimates 

In the next step we have to estimate the eigenvalues of L^ and — A^. Let us 
start with the longitudinal part. It is easy to check the identity 

L± = Ct{d)S±(c v + C h ^jd- 

combining it with the minimax principle and the fact that the eigenvalues of 
S behave as (^-) 2 £ 2 + 0(1) as I — > ±00, we arrive at the following conclusion: 



Lemma 3.5 There is a positive C such that the eigenvalues lf(d) of Lj, 
numbered in the ascending order, satisfy the inequalities 

\lf(d) -nd < Cj 2 d (3.5) 

for all j G N and d small enough. 

The transverse part is a bit more involved. Our aim is to show that in 
the strong-coupling case the influence of the boundary conditions is weak, 
i.e. that the negative eigenvalues of the operators — A^ do not differ much 
from the number (|2.7j) . 

Lemma 3.6 There exist positive numbers Ci, 1 < i < 4, such that each one 
of the operators — A^ has exactly one negative eigenvalue t„ which satisfies 

£ Q - S(a) <t-<Z a <ti<£ a + S(a) (3.6) 

for a large enough negative, where 

S(a) := C 1 C 2 y/dC e exTp(-C 2 dCa) 
with ( a := (— ^a) 1 / 2 , provided d( a > C3 and dM < C4. 
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Proof. Let us start with the eigenvalue of the operator — A+ involved in the 
upper bound; the argument will be divided into four parts. 

1. step: We will show that the number — k 2 a with k a > is an eigenvalue 
of — A+ iff k a is a solution of the equation 

x = ( a r](x) , (3.7) 

where ( a has been defined above and r\ is the function given by 

„: R+ ^ H+ , , W = exp(-|M); (3.8) 

the symbols K Q , J denote the Macdonald and the modified Bessel function, 
respectively |ASj . To verify this claim we note that the eigenfunction ip of 
— A+ corresponding to — k 2 a is a linear combination 

<p(r) = DMk a ry/ 2 + D 2 K,{k a r)r 1 ' 2 

with the coefficients D\, D 2 chosen in such a way that the conditions follow- 
ing from Lp~ D.bc(dBd) and r~ 1 ' 2 ip~a.bc(0) are satisfied. Using the behavior 
of K , I at the origin 

K (p) = -\n^ + ^(l) + O(p) and J Q (p) = 1 + 0{p) , (3.9) 

as p — > 0, we can readily check that ip fulfils the needed boundary conditions 
iff (Di, D 2 ) G ker M(a), where M(a) is the matrix given by 

m,» = ( h{ Y ] K ? {k t d \ ) 

3K ' V 1 uj{a,k a ) J 

with u>(a, k a ) := tp(l)—27ia—hi(k a /2). Of course, the condition kerM(a) ^ 
is equivalent to det M(a) = 0; the latter holds iff k a is a solution of ()3.7|) . 

2. step: Our next aim is show that the equation (j3.7J) has at least one 
solution for —a sufficiently large, and moreover, that such a solution k a 
satisfies the inequalities 

C(a <K<Ca (3.10) 

with C G (0, 1) independent of a. Using again ()3.9|) together with the asymp- 
totic behavior of the functions K , 1 at infinity, we get for a fixed a 

(aVi X ) ~ y (a as X — » OO 
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and 

CaV( x ) = 9a,d% + 0(x 2 ) as x^O, (3-11) 

where g a> d '■= ^e~^^d( a . It is clear that the error term is uniform with 
respect to a over finite intervals only, however, if 

9a,d > 1 (3.12) 

then the equation (|3.7|) has obviously at least one solution. The second 
inequality in (|3.1U|) holds trivially because T](x) < 1 for any x > 0. Let us 
assume that the first one is violated. This means that there is a sequence 
{a n } with a n — > — oo as n — > oo such that i](k an ) — > as n — > oo. This 
may happen only if the fc an tends to the singularity of K , in other words if 
k an — > holds as n — > oo. However, the inequality (|3.12|) is valid for a n with 
n large enough, thus small k an can not in view of the asymptotics (j3.11|) be 
a solution of (|3.7J) in contradiction with the assumption. 

3. step: To show that there exists only one solution of ()3.7)1 it suffices to 
check that the function h a : ¥L + i— > R, 

= x - Ca^(^) , 

is strictly monotonous for x G (C'Ca; (a) and —a sufficiently large. Using 
again the behavior of Kq, Iq at large values of the argument we find that the 
derivative rf(x) —> as x — > oo which implies the result. 

4- step: It remains to show that the eigenvalue t+ = — k 2 a satisfies the 
second one of the inequalities 

ta< -hi <U + S{a). (3.13) 

Since the functions —K , I are increasing and io(0) = 1 we get from (|3.10j) 
the estimate 

V(K) > exp (-ifo(C'Cad)) • 

Putting now S(a) — (l — exp ^— 2K (C( a d)j^j ( 2 and using the asymptotic 
behavior of Kq at large distances one finds that 

S(a) < C'lCaV^Ca ex P( — C*2^Ca) as a ~" ^ — 00 

holds with suitable constants C\, C 2 and the inequality ()3.13|) is satisfied 
which concludes the proof for the operator — A+. 
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Let us turn to the operator — A~. The argument is similar, so we just 
sketch it with the emphasis on the differences. The number t~ = —k 2 a is an 
eigenvalue of — A~ iff k a is a solution of the equation 

x = ( a fj(x), (3.14) 

where fj : M + — > M + is the function given by 

(S (xd?\ \ 
— tt~, — rr ) 5 S F (xd) = FAxd^xd + w d F (xd) 
Si{xd) J 

for F = K, I, where I\ = Ii, K\ — —K\ and Wd ■= \ — Md\ we assume that 

w d >0. (3.15) 

To proceed further, we employ again the asymptotics of functions /„, K n , 
n = 0, 1, for x — > and at large values of the argument. It is easy to see 
that the behavior of x t— > for small x is dominated by that of -Ko(')- 

Thus mimicking the second step of the above argument we can show that the 
equation ()3.14|) has at least one solution for —a sufficiently large provided 
that assumption (|3.15|) is satisfied. Repeating the third step we can check 
that the solution k a is unique for —a sufficiently large. By reduction ad 
absurdum, as in the second step, we can also prove that there exists C such 
that C( a < k a , which means that k a — > oo as a — >• — oo. The constant can be 
made more specific: using the fact that the term — dxK\(dx) dominates the 
behavior of Sx(dx) for large x and Si > we get fj > 1 for —a sufficiently 
large, i.e. 

Using properties of the special functions involved here we also find that 

fj(k a ) < exp(CKi(dC a )dC a ) 

holds for any C satisfying {wd)~ l + (d(a)^ 1 < C . Thus proceeding similarly 
as in the fourth step we infer that there are constants C±, C2 such that 

£a (3-16) 

where 

S(a) < CiCq \fd~Ca ex p( - C*2^Ca) as a — * —°° ■ 
Finally putting together (J3~T2|) . (j3~T3|l . (j3T5) and (jHHEj) we get the claim 
with G\ := max{Ci,(7i} and Ci := min{C2,C72}. ■ 
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3.4 Proof of Theorem 13.31 for a loop 

Suppose now that T is a closed curve. The result will follow from combination 
of the above estimates. We have to couple the width of the neighbourhood 
fid and the coupling constant a in such a way that d shrinks properly to zero 
as a — > — oo. This is achieved, e.g., by choosing 

d(a)=e™. (3.17) 

Proof of (a): To find the asymptotic behavior of eigenvalues Xj(a) of H a ^ 
we will rely on the decomposition (|3.4j) . according to which we know that 
the negative eigenvalues of H a r are squeezed between lf(d) +t„. Since the 
operators — have a single negative eigenvalue, the sought values Xj(a) are 
ordered in the same way as lf(d) are. Combining ()3.17|) with the results of 
Lemmas 13.51 and 13.61 we get for the upper and lower bound 

lf(d(a)) +t± = £ a + ft + 0{f a ) as a -> -oo , 

and of course, the same asymptotics holds for Xj(a). Clearly, to a given 
integer n there exists a(n) el such that l+(d(a)) + t+ < is true for all 
a < a(n); this completes the proof of (a). 

Proof of (b): Using the above asymptotic estimates and Lemma f3. 51 we get 

vj{a) < Xj{a) < u+(a) , (3.18) 

where 

vf(a) := ^ + f ( ( 2 + ( e ™)) ± v 

and v = 4~ 1 ||k 2 || 00 . Combining this with the minimax principle we arrive at 
the two-sided estimate 

${j e Z : u+(a) < 0} < ta d (H a ) < ${j e Z : vj{a) < 0} , 

which implies 

M^«) = -(-e«») 1/2 (n-o(e TO )). 

7T 
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3.5 A curve with free ends 

The part (b) of Theorem 13 . 31 does not require T to be a closed curve. One can 
repeat the argument with a small modification taking for Qd a closed tube 
around T bordered by the additional "lid" surfaces normal to T at its ends. 

j1 2 

Thus instead of S we have a pair of comparison operators S l = — — ^ on 
L 2 (0, L) with D{S i ) := {/ G W 2 ' 2 (0, L), f~i.bc}, i — D, N, which give in the 
same way as above estimates for the eigenvalues Aj(a) of H a r as a — > — oo, 
namely 

£ Q + /if + 0(O < A» <£ a + + 0(e™), 

where /i*, j = 1, 2, . . . , denote the eigenvalues of 5**. The fact that the latter 
are different for the Dirichlet and Neumann condition does not allow us to 
squeeze Xj(a) sufficiently well to get its asymptotics in analogy with the 
claim (a) of the theorem. On the other hand, the behavior of yJ^ as 
j — >• oo allows us to find an asymptotic estimate for the counting function. 
Recall that the eigenvalues of —A' = — : D(S l ) — > L 2 (0,L) are of the 
form s* = j 2 {j) 2 , where j G N for i = D and j G N U {0} for i — N; thus in 
analogy with (|3.18|) we can define the functions 

4(a) :=a+jl((^) 2 + 0(e 7ra ))±^ 

where j + = j and j~ = j — 1 with j G N, which give a two-sided bound 
for Xj(a). Combining it with the minimax principle we arrive again at the 
formula 

M# a , r ) = -(-U 1/2 (l + 0(e™)) as «--oo. 



Remark 3.7 While Theorem 13.31 was formulated for a single finite curve, 
which may not be closed for part (b), the argument easily extends to any 
T which decomposes into a finite disjoint union of such curves, up to the 
eigenvalue numbering. The latter may be ambiguous in case that the corre- 
sponding operator S, which is now an orthogonal sum of components of the 
type (|3.ip . exhibits an accidental degeneracy in its spectrum. 
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4 Infinite asymptotically straight curves 

We know from jEKj that the operator H a Y has a nonempty discrete spectrum 
if T is an infinite C 4 curve which is non straight but it is asymptotically 
straight in the following sense 

(ar inf l) for all s G R we have \k(s)\ < M|s| _/3 , where (3 > 5/4 and M > 0. 
Moreover one has to assume that 

(ar in f2) there exists a constant c G (0, 1) such that \j(s) —j(s') | > c\s — s'\. 
If these conditions are satisfied then the operator H a p is self-adjoint and 

(Ha,r) = [£<*, oo) , cr d (iJ Qj r) 7^ . 

Since the infinite curve has no free ends, the asymptotics of eigenvalues of 
H a T for a — ► — oo can be found in the same way as for the loop. We employ 
the comparison operator which now takes the form 

S = ~ ~ \<sf : D(S) - L 2 (R) 

with the domain D(S) equal to W 2 ' 2 (R). It is a Schrodinger operator on line 
with a potential which is purely attractive provided k^O, and therefore 

a d (5)^0. 

On the other hand, in view of the assumed decay of curvature as \s\ — > oo 
the number N := ^a^(S) is finite |RS| Thm. XIII. 9]. Using the symbol fij 
for the j-th eigenvalue of the operator S we get the following result. 

Theorem 4.1 Under the above stated assumptions there is «o G R such that 
tjo"d(-^a,r) — N holds for all a < «o- Moreover, the j-th eigenvalue Xj(a) of 
H a ,r, j — 1, • • • j N, admits the asymptotic expansion 

Xj(a) = £ a + //j + (^(e™) as a — > — oo . 
Since the proof fully analogous to that of Theorem 13.31 we omit details. 
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5 Spectrum for an infinite periodic curve 

5.1 The Floquet— Bloch decomposition 

Now we turn our attention to Hamiltonians with singular perturbations sup- 
ported by a periodic C 4 curve without self-intersections. In other words we 
assume that there is a vector K x = K G R 3 and a number L > such that 



Of course, we can always choose the Cartesian system of coordinates such 
that K = (K, 0, 0) with K > 0, and 7(0) = 0. As usual in periodic situations 
we decompose the space R 3 according to the periodicity of T. To this aim 
we define the basic period cell as 



where {Kj} 3 =1 are linearly independent vectors in R 3 ; without loss of general- 
ity we may suppose that K 2 _L K 3 . Then the translated cells C n := C + nK, 
where n G Z, are mutually disjoint for different values of the index and 
R 3 = [j neZ C n . As in the previous section we assume that T has no self- 
intersections. However, to proceed further we need an additional assumption, 
namely 

(ar per ) the restriction of Tq := C fl T to the interior of C is connected. 

Let us note the choice of the point s = is important in checking the as- 
sumption (ar per ), and for the same reason we do not require generally that 
Ki _L {K 2 ,K 3 } (see also Remark 15.41 below) . 

While a smooth periodic curve without self-intersections satisfies (an), 
the property (ar per ) ensures that we can choose a neighbourhood of T c which 
is connected set contained in C; this is important for the construction de- 
scribed below. In view of Theorem 12.31 the Hamiltonian with the singular 
perturbation supported by T is well defined as a self-adjoint operator in 
L 2 (R 3 ). To perform the Floquet-Bloch reduction for H a Y we decompose 
first the state Hilbert space into a direct integral 



7(5 + L) = K + 7 (s) for all s G R . 



3 




(5.1) 



8=1 
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It is a standard matter to check that the operator U : L 2 (R 3 ) — > H given by 
(Uf)e(x) = -4^ E e^VO* + nK) (5.2) 

on / G C£°(R 3 ) acts isometrically, so it can be uniquely extended to a unitary 
operator on the whole L 2 (M 3 ). We will say that the function / G C 2 (C\rc) 
belongs to T a (0) if it satisfies the condition 

f~a.bc(T e ) , 

and furthermore, for all x such that both x and x + K belong to dC and 
x 7^ (0, 0, 0) we have 

fM(x + K) = e» K fM(x), i/ = 0,l, (5.3) 

where J*- -* := /, := 9 X1 /. Now we define H a ^{9) as the self-adjoint 
Laplace operator in L 2 {C) with the boundary conditions introduced above; 
more precisely, H a (9) is the closure of 

H a ,r(9) : D(H a , r (6)) = {/ G T Q (#) : H a ^(9)f G L 2 (C)} — > L 2 (C) , 
H a , T (e)f(x) = -Af(x), xeC\T c . 

The following lemma states the usual unitary equivalence between if Qi r and 
the direct integral of its fiber components H a> r(9). 

Lemma 5.1 UH^U' 1 = J { % K ^ /K) H a , r {6) d9. 

Proof. Take a function / belonging to the set 

C ■= {g G C 2 (M 3 \r) : f~a.bc(T) , supp / is compact } (5.4) 

then for alH = 1, 2, 3 we have 

(Ud i f(x)) = d i (Uf) (x), x£T, 

and the same relations hold for the second derivatives. Thus to prove the 
lemma it suffices to show that any function admitting the representation 
(U f)e with / G C belongs to T a (9). It is easy to check that for all x ^ (0, 0, 0) 
such that x and x + K are in dC we have 

{(Uf)$\x + K) = e» K ((Uf)$\x)) for v = 0, 1 . 
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The behavior of the function {Uf)g in the vicinity of Tq is characterized by 
the limits E((Uf)g)(-) and Q((U f)g)(-). Using the periodicity of T we get 

Z((Uf)e)(s) = (2 7 r)- 1 /^e-^H(/)( S + nL), s G (0, L) , 
mUf)e){s) = (2 7 r)- 1 /^ e — "<a(f)(s + nL) , s G (0, L) ; 

to derive these relations we used also the uniform convergence of the sums. 
In this way we conclude that (Uf)g~a.bc(Tc). The Laplace operator in L 2 (C) 
with the domain consisting of functions which admit the representation (Uf)$ 
with / G C is essentially self-adjoint and its closure coincides with H a ^(0); 
this completes the proof. ■ 

5.2 Spectral analysis of H a ^(0) 

As in the case of a finite curve we can now analyze the discrete spectrum of 
the operator H a ^(0). Before doing that let us localize the essential spectrum. 
An argument analogous to that of Proposition 13.21 shows that the singular 
perturbation supported by Tc does not change the essential spectrum of the 
Laplacian in a slab with Floquet boundary conditions, i.e. 

a ess (H a , r (6)) = [6\oo) . (5.5) 

To describe the asymptotic behavior of the eigenvalues of H a ^{6) we in- 
troduce a comparison operator by Sg = — 4^ — ^j- : D(Sg) — > L 2 (0,L), 
where 

D(S e ) := { / G W 2 > 2 (0, L) : f(L) = e i6K f(0), f{L) = e ieK f'(0) } . 
In analogy with Theorem 13.31 we state: 

Theorem 5.2 Under the assumption given above for a fixed number n there 
exists a(n) G R such that §ad(H a (6)) > n holds for a < a(n). Moreover, the 
j-th eigenvalue of H a ^{6) has the asymptotic expansion of the form 

Xjia^) =t a + iij{0) +0{e™) as a^-oo, 

where fij{0) is the j-th eigenvalue of Sg and the error term is uniform with 
respect to 9. 
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Proof. The argument follows closely that of Theorem 13 .3| the only difference 
is the replacement of periodic boundary condition by the Floquet one. The 
fact that the error is uniform w.r.t. 9 is a consequence of Lemma 13.51 and 
continuity of the functions /^(-). ■ 

5.3 Spectral analysis of H a ^ in terms of H a ^{9) 

Now our aim is to express the spectrum of H a p in the terms of H a ^{9). First, 
let us note that combining (|5.5|) with standard results |RS| Sec. XIII. 16] we 
get the following equivalence for the positive part of spectrum 

a{H a , T ) n [0, oo) = |J a(F a , r (0))n[O,oo) = [O,oo). 

ee[-7r/K,n/K) 

The negative part of spectrum is more interesting being given by the union of 
ranges of the functions Xj{a, •). They give rise to well-defined spectral bands 
because the latter are continuous in the Brillouin zone [—7t/K,tt/K). This 
can be seen by checking in the usual way, putting 9 into the operator and 
showing that the 9 dependent part is an analytic perturbation. Alternatively, 
one can take g = (Uf)$ with / e C as defined by (|5.4jl and investigate the 
functions 

qg{0) := (g,H a , r (e)g) L 2 {c) = i- ]T e-« n -^ 9 (f n , H a , r f m ) LHc) , 

n,m£Z 

where f n {x) '■— f{x + riK). In view of (|5.2j) and the uniform convergence of 
the respective sums such a q g (-) is continuous for g runing over a common 
core of all H a ^{9). Thus by the minimax priciple we get the continuity of 
Xj(a, •) and combining this fact with the results of |RSj we get 

a(H a>r ) n (-oo, 0] = I |J o{H a {9)) J n (-oo, 0] 

arriving finnaly at 

tr(F a) r) = |J a(H a>r (9)) . 

0e[-Tr/K,n/K) 

These results together with Theorem 15 . 21 allow to describe the band structure 
of H a< r, in particular, the existence of gaps. Notice that this operator as well 
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as S = -4 - in £ 2 (IR) commute with the complex conjugation, so their 

Floquet eigenvalues are generically twice degenerate depending on \6\ only. 
For the comparison operator thus width of the j— th gap is 

G,(S) = (i j+1 (n/K) - Hjfr/K) for odd j 
/"i+i( ) - for even 3 

and similarly for H a Y- The expansion of Theorem 15.21 then gives 

G 3 (H a ^) = G 3 (S) + 0(en- 

In combination with the known result about existence of gaps for one-dimen- 
sional Schrodinger operators we arrive at the following conclusion. 

Corollary 5.3 Suppose that in addition to the above assumption the function 
k(-) is nonconstant. In the generical case when S has infinitely many open 
gaps, one can find to any n G N an a(n) G M such that the operator H a Y 
has at least n open gaps in its spectrum if a < a(n). If the number of gaps 
in o~(S) is N < 00, then o~(H a Y) has the same property for —a large enough. 

Notice that this property is determined by the curvature alone. Thus the 
result does not apply not only to the trivial case of a straight line, but also 
to screw-shaped spirals V for which k is nonzero but constant. 

Remark 5.4 It is not always possible to choose C in the form of a rectan- 
gular slab (|5.1|) as we did above, which would satisfy the assumption (ar por ); 
counterexamples can be easily found. However, if we choose instead another 
period cell C with a smooth boundary for which the property (ar per ) is valid, 
the argument modifies easily and the claim of Theorem 15.21 remains valid. 
On the other hand, such a decomposition may not exist if the topology of 
T is non-trivial; a simple counterexample is given by a "crotchet- shaped" 
curve. While we conjecture that the claim of Theorem 15.21 is still true in this 
situation, a different method is required to demonstrate it. 

5.4 Compactly disconnected periodic curves 

So far we have considered a single periodic connected curve. A slightly 
stronger result about the existence of gaps in spectrum of H a r as a — » —00 
can be obtained for compactly disconnected periodic curves in M 3 , i.e. such 
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that they decompose into a disjoint union in which each of the connected 
components is compact. To be more specific, we consider a family of curves 
obtained by translations of a loop Tq (being a graph of a function 70) gener- 
ated by an r-tuple {Kj} linearly independent vectors, where r = 1, 2, 3. The 
curve T in question is then a union T = IJnez r ^n, where T n are graphs of 

7n := 7o + HiKi : [°> L ] ~" R3 ' n = ; 

for the sake of brevity we put here T no = r , 7 no = 70, where no := (0, 0, 0). 
We assume that T is contained in the interior of the period cell 

C= |X> K * : <U < lj x {K;}\ 

which is noncompact if r = 1,2 and compact otherwise. Similarly as before 
we can make Floquet-Bloch decomposition of H a T into a direct integral of 
the fiber operators H a ^(9). However now, since dist(<9C,r c ) > holds by 
assumption, the comparison operator S = S(0) is now independent of the 
quasimomentum 9 G ]^[^ < < ^ <r .[ — TrlK^I" 1 , ttIK^I" 1 ). While in the previous case 
some gaps of S(9) might be closed, now they are all open. As a result each 
gap in the spectrum of a(H a r ), which depends of course on 9, will eventually 
open for —a large enough. 

Theorem 5.5 Under the assumptions stated above the spectrum of H a j<(9) 
is purely discrete if r = 3, and a css (H a ^(9)) = [X^[=i @h 00 ) if r = 1> 2. The 
j-th eigenvalue of H a ^{9) admits the asymptotic expansion of the following 
form, 

\j(a, 9) = £ a + /ij + 0(e Ka ) as a — > —00 , 

where fij is the j-th eigenvalue of S and the error is uniform w.r.t. 9. Con- 
sequently, for any n G N there is a(n) G M such that the operator H a ^ has 
at least n open gaps in its spectrum if a < a(n). 

6 Concluding remarks 

(a) The results obtained in the previous discussion can be rephrased as a 
semiclassical approximation. To see this let us consider the Hamiltonian 
H a ,r(h) with the Planck's constant h reintroduced; the latter is understood 
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in the mathematical sense, i.e. as a parameter which allows us to investigate 
the asymptotic behavior as h — * 0. The operator in question then acts as 

H a , r (h)f(x) = -h 2 Af(x) , x € R 3 \ T , 

and has the domain 

D{H a , T {h)) = {f G T M a : f~a(h).bc(r)} , 

where ^ 

a(h) := a H ln/i. (6.6) 

27T 

This definition of H a ^{h) requires a comment. In the case codimT = 1 
discussed in |Exj the Hamiltonian is defined by the natural quadratic form, 
hence introducing h means a multiplicative change of the coupling parameter, 
a — > ah~ 2 ; one can see that also from the approximation of such an operator 
by means of scaled regular potentials [Elj . 

In contrast to that a two-dimensional point interaction involves a com- 
plicated nonlinear coupling constant renormalization AG HH[ Sec. 1.5], so 
introducing Planck's constant is in this case arbitrary to a certain extent. 
We choose the simplest way noticing that the relation between the free oper- 
ators —A and — h 2 A can be expressed by means of the scaling transformation 
x i — > hx, and require the similar behavior for the singular interaction term; 
it is well known that a scaling for a two-dimensional point interaction is 



equivalent to a logarithmic shift of the coupling parameter - cf. EG STj. In 
view of (|6.6|) the semiclassical limit h — > is within this convention for a 
fixed coupling constant a equivalent to a(h) — > — oo which means a strong 
coupling again. Since H a ^[h) = h 2 H a ( h ^ r (l) we see that the eigenvalues 
Aj(a, h) of H a ^(h) take then the following form, 

Aj(ar, h) = £ a + fijh 2 + 0(h 5/2 ) as h -»• . 

In the same way we find the counting function which is given by 

#a d (H a (h)) = 4(-U 1/2 a + 0{h^)) . 
Tin 

(b) Let us finally list some open problems related to the present subject: 
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One is naturally interested in the asymptotic expansion in the situation 
when T is a curve with free ends and the present method allows us to 
treat the counting function only; the analogous question stands for 
planar curves |EYlj and surfaces with a boundary |Exj . We conjecture 



that the expansion of Theorem 13.31 holds again with \ij corresponding 
to the comparison operator which acts according to (J3.1|) with Dirichlet 
boundary conditions at the boundary of T. 

The results can be extended to higher dimensions provided codim T < 3 
so that the singular interaction Hamiltonian is well defined. 

The smoothness assumption is crucial in our argument. A self-similar 
curve such as a broken line consisting of two halflines joined at a point 
provides an example of a situation where the asymptotic behavior dif- 
fers from that of Theorem 13.31 One can ask, e.g., how the asymptotics 
looks like for a piecewise smooth curve with non-zero angles at a dis- 
crete set of points. 

Another important question concerns the absolute continuity of the 
spectrum in case when T is a periodic curve or a family of curves. The 
answer is known if codim T = 1 and the elementary cell is compact 
. The cases of a single connected periodic curve or a periodic 



surface diffeomorphic to the plane are open, and the same is true for 
periodic curve(s) in M 3 , i.e. the situation with codim T = 2. 
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